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Abstract 

The relativistic mean-field theory provides a framework in which the nuclear many- 
body problem is described as a self-consistent system of nucleons and mesons. In 
' the mean-field approximation, the self-consistent time evolution of the nuclear system 

. describes the dynamics of collective motion: double giant resonances, nuclear com- 

I pressibility from monopole resonances, regular and chaotic dynamics of isoscalar and 

' isovector collective vibrations. 

1 Introduction and outline of the model 

^ I Relativistic mean-field (RMF) models have been successfully applied in calculations of nu- 
clear matter and properties of finite nuclei throughout the periodic table. In the self- 



consistent mean-field approximation, detailed calculations have been performed for a variety 
of nuclear structure phenomena In the present work we review the applications of RMF 
', to the dynamics of collective vibrations in spherical nuclei. In relativistic quantum hadro- 
dynamics , the nucleus is described as a system of Dirac nucleons which interact through 
the exchange of virtual mesons and photons. The Lagrangian density of the model is 
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- gaipaip - g^i/j'j ■ uijj - Qp'ip'^ ■ prip - e^7 ■ A^^^—^^ip . (1) 

The Dirac spinor denotes the nucleon with mass m. ma, m^, and uip are the masses of the 
cr-meson, the w-meson, and the p-meson, and Qa, g^j, and gp are the corresponding coupling 
constants for the mesons to the nucleon. U{a) denotes the nonlinear a self-interaction 

U{cr) = \mla^ + \g,a^ + \g,a\ (2) 
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and VL^''', R^"', and F^"' are field tensors ||. 

The coupled equations of motion are derived from the Lagrangian density (|1|). The Dirac 
equation for the nucleons: 



cx (-iV - Qi^uj - gpTp - e-^^— J + /9(m + g^a 



+gcoUJo + gpTpo + e Ao 



2 

(3) 



and the Klein-Gordon equations for the mesons: 

[d^-A + ml)a = -g^ps- g2(T^ - g-ifT^^ (4) 

(d^ - A + ml)ujp = gjp (5) 

{d^-A + ml)p, = gj, (6) 



5* -AM. = ejT. (7) 



In the relativistic mean-field approximation, the nucleons described by single-particle spinors 
tpi {i = 1,2,..., A) are assumed to form the A-particle Slater determinant |$), and to 
move independently in the classical meson fields. The sources of the fields, i.e. den- 
sities and currents, are calculated in the no-sea approximation [Q: the scalar density: 
Ps = Z^iLiV'iV'i) the isoscalar baryon current: = J2t=i''Pi'y^ipi, the isovector 

baryon current: = J2t=i4'i1^^''Pi^ the electromagnetic current for the photon-field: 

Jem = ^t=i i'll^^-^i'i- The summation is over all occupied states in the Slater determi- 
nant 1$). Negative-energy states do not contribute to the densities in the no-sea approxi- 
mation of the stationary solutions. It is assumed that nucleon single-particle states do not 
mix isospin. 

The ground state of a nucleus is described by the stationary self-consistent solution of the 
coupled system of equations (H)-®, for a given number of nucleons and a set of coupling 
constants and masses. The solution for the ground state specifies part of the initial conditions 
for the time-dependent problem. The dynamics of nuclear collective motion is analyzed in 
the framework of time-dependent relativistic mean-field model, which represents a relativis- 
tic generalization of the time-dependent Hartree-Fock approach. For a given set of initial 
conditions, i.e. initial values for the densities and currents, nuclear dynamics is described by 
the simultaneous evolution of A single- particle Dirac spinors in the time- dependent mean 
fields. Frequencies of eigenmodes are determined from a Fourier analysis of dynamical quan- 
tities. In this microscopic model, self-consistent time-dependent mean-field calculations are 
performed for multipole excitations. An advantage of the time-dependent approach is that 
no assumption about the nature of a particular mode of vibrations has to be made. Retar- 
dation effects for the meson fields are not included in the model, i.e. the time derivatives df 
in the equations of motions for the meson fields are neglected. This is justified by the large 
masses in the meson propagators causing a short range of the corresponding meson exchange 
forces. Negative energy contributions are included implicitly in the time-dependent calcu- 
lation, since the Dirac equation is solved at each step in time for a different basis set 0. 
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Negative energy components with respect to the original ground-state basis are taken into 
account automatically, even if at each time step the no-sea approximation is applied. 
The description of nuclear dynamics as a time-dependent initial-value problem is intrinsi- 
cally semi-classical, since there is no systematic procedure to derive the initial conditions 
that characterize the motion of a specific mode of the nuclear system. The theory can be 
quantized by the requirement that there exist time-periodic solutions of the equations of 
motion, which give integer multiples of Planck's constant for the classical action along one 
period ||^. For giant resonances the time-dependence of collective dynamical quantities is 
actually not periodic, since generally giant resonances are not stationary states of the mean- 
field Hamiltonian. The coupling of the mean-field to the particle continuum allows for the 
decay of giant resonances by direct escape of particles. In the limit of small amplitude os- 
cillations, however, the energy obtained from the frequency of the oscillation coincides with 
the excitation energy of the collective state. In Refs. 0, ^, ^ we have shown that the model 
reproduces experimental data on giant resonances in spherical nuclei. 

2 Dynamics of collective vibrations 
2.1 Double Giant Resonances 

The physics of giant resonances in nuclei has been the subject of extensive experimental 
and theoretical studies for many years. However, only recently substantial evidence has 
been reported for two-phonon states built with giant resonances. Resonant structures that 
were observed in heavy-ion inelastic scattering, have been interpreted as possible multiple 
excitations of the giant quadrupole resonance. Double giant dipole resonances have also been 
discovered in the neutron and in the 7-spectra of nuclei that have been Coulomb excited in 
relativistic heavy-ion collisions. In Ref. ^ we have performed time-dependent relativistic 
mean-field calculations and found evidence for modes which can be interpreted as double 
resonances, and which in a quantized theory correspond to two-phonon states. 
As an example of double resonances in light nuclei, we consider the double isoscalar giant 
quadrupole resonance in ^''Ca. The experimental spectrum exhibits a prominent structure, 
centered at 34 ± 2 MeV excitation energy, with a width of 9 ± 2 MeV. It is interpreted as 
the two-phonon state of the single isoscalar GQR at 17.5 MeV. Using the NL-SH parameter 
set for the effective Lagrangian, we have studied isoscalar quadrupole oscillations in ^°Ca. 
The quadrupole mode of oscillations is excited by deforming the spherical solution for the 
ground-state. For a specific initial deformation, we follow the time-evolution of the collective 
variable, the quadrupole moment 

q2o{t) = {m\Q2o\m) = {m\r'Y2omt)) (s) 

The time-dependent quadrupole moment shown in Fig. 1 corresponds to an initial axial 
deformation of the baryon density (3 = 0.38. The resulting Fourier spectrum displays a 
strong peak at 18.5 MeV, in reasonable agreement with the experimental data. We find 
evidence for excitation of a higher mode in the oscillations of the baryon density. The wave 
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function of the nuclear system is a Slater determinant at all times, and therefore can be 
expanded in the basis of the ground state |$o) 



<^{t)) = \%) +J2zmi{t)a+ai\%) + z^i{t)z„,>iit)a^a+,ai,ai\%) + ... (9) 



If the total wave function contains collective 2p — 2h components, they will be observed in 
the Fourier spectrum of the time- dependent baryon density pelr, t) = Y.t=i t)ipi{r, t). 
Because of axial symmetry and the isoscalar nature of the excitation, it is sufficient to 
consider oscillations in time of the baryon density on the positive z-axis. In Fig. 2 we 
display the Fourier transforms of the time- dependent baryon density for various values of 
the coordinate z. Two peaks are clearly observed. The first one at 18.5 MeV corresponds to 
the isoscalar quadrupole resonance. It gradually increases from the center toward the surface 
of the nucleus. If we plot the values of the Fourier transforms at 18.5 MeV as a function of z, 
the resulting curve corresponds to the transition density. The transition density for the first 
peak is typical for isoscalar quadrupole resonances. The second peak is at 37 MeV, twice 
the energy of the GQR. It has a maximum in the center of the nucleus, at first decreases 
with z, but then appears again on the surface. Compared to the GQR transition density, 
the curve for the 37 MeV peak displays an additional node. 

Periodic solutions of the time-dependent Dirac equations can be used to construct the ener- 
gies of the many-body system. The energy spectrum can be obtained from a semi-classical 
quantization procedure. One finds periodic solutions such that the mean-field action along 
a periodic orbit 



is equal to an integer multiple of the Planck constant / = nh. In Eq. (0) T is the period 
of oscillations, > denotes time-dependent single- nucleon Dirac spinors, are the 

corresponding single-nucleon energies in the unperturbed ground-state, and the summation 
runs over the occupied states. Because of the coupling to the continuum in the mean-field 
description, giant resonances are not stationary states of the Hamiltonian. Consequently, 
a non-periodic dependence on time is obtained for dynamical quantities. If the damping is 
very strong, the giant resonance is not periodic even on the average, and the quantization 
condition cannot be applied. However, if the motion is nearly periodic, i.e. the damping is 
relatively weak, the quantization procedure can still be used to calculate the energies, and 
the effect of damping can be taken into account approximately. 

In Fig. 3 we display the mean-field action as function of the excitation energy of the nucleons 
and of the initial deformation /?. For the values of (3 indicated by dots we have integrated 
the coupled system of Dirac and Klein-Gordon equations and calculated the action integral. 
The mean-field action is a quadratic function of the initial deformation P, and an almost 
perfect linear function of the excitation energy. Only above 60 MeV a slight deviation from 
a pure linear dependence is observed. The values of < -E* > for which the action is an 
integer multiple of the Planck constant are: I = 1 h for 18.5 MeV, and I = 2 h for 37.1 
MeV. A double giant dipole resonance has been observed in relativistic Coulomb excitation 
of 208pb. The single GDR is found at 13.3 ± 0.1 MeV with a width of 4.1 ± 0.1 MeV. 
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The sum energy of coincident photon pairs displays a broad structure at 25.6 ± 0.9 MeV 
with a width of 5.8 ± 1.1 MeV. It is interpreted as the double GDR. In order to excite 
isovector dipole motion we define the initial conditions: at t = (in the center of mass 
system) all protons start moving in the +z direction with velocity v.,^, and all neutrons 
start moving in the —z direction with velocity Vu = j^v-j^. For the NL-SH parameter set, 
the Fourier spectrum of the time-dependent dipole moment displays a strong peak at 12.9 
MeV excitation energy, in good agreement with experimental data. We have found that the 
mean-field action is an integer multiple of the Planck constant: I = 1 h for < E* >= 12.9 
MeV, and I = 2 h for < E* >= 25.9 MeV. Therefore, the energy of the one-phonon state, 
calculated from the mean-field action, coincides with the resonant energy of the mean peak 
in the Fourier spectrum, and the two-phonon state at 25.9 MeV is in excellent agreement 
with the experimental value for the excitation energy of the double GDR. 

2.2 Monopole Giant Resonances and nuclear compressibility 

The study of isoscalar monopole resonances in nuclei provides an important source of in- 
formation on the nuclear matter compression modulus i^nm- This quantity is crucial in 
the description of properties of nuclei, supernovae explosions, neutron stars, and heavy-ion 
collisions. In principle the value of can be extracted from experimental energies of 
isoscalar monopole vibrations in nuclei (giant monopole resonances OMR). However, the 
complete experimental data set on isoscalar GMR does not limit the range of K^^^ to better 
than 200 — 350 MeV. Microscopic calculations of GMR excitation energies might provide 
a more reliable approach to the determination of the nuclear matter compression modulus. 
Modern non-relativistic Hartree-Fock plus RPA calculations, using both Skyrme and Gogny 
effective interactions, indicate that the value of should be in the range 210-220 MeV. 
In relativistic calculations on the other hand, both time-dependent and constrained RMF 
results indicate that empirical GMR energies are best reproduced by an effective force with 
i^nm ~ 250 - 270 MeV. 

In Ref. we have performed time-dependent and constrained RMF calculations for 
monopole giant resonances of a number of spherical closed shell nuclei, from ^^O to the 
heavy nucleus ^°^Pb. For the effective Lagrangian we have used six parameter sets, which 
differ mostly by their prediction for Knm, but otherwise reproduce reasonably well experi- 
mental data on nuclear properties. The idea is to restrict the possible values of the nuclear 
matter compression modulus, on the basis of the excitation energies of giant monopole states 
calculated with different effective interactions. In addition to the four non-linear sets NLl, 
NL3, NL-SH and NL2, we have also included two older linear parametrizations, HS and LI. 
The sets NLl, NL-SH and NL2 have been extensively used in the description of properties 
of finite nuclei [|I|. In order to bridge the gap between NLl {K^m = 211.7 MeV), and NL-SH 
{Knm = 355.0 MeV), we have also included a new effective interaction NL3 |^ {K^m = 271.8 
MeV). This new parameter set provides an excellent description not only for the properties 
of stable nuclei, but also for those far from the valley of beta stability. From the energy 
spectra and transition densities calculated with these effective forces, it has been possible 
to study the connection between the incompressibility of nuclear matter and the breathing 
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mode energy of spherical nuclei. For the isoscalar mode we have found an almost linear 
relation between the excitation energy of the monopole resonance and the nuclear matter 
compression moduhis. For the determination of especially relevant are microscopic 
calculations of GMR excitation energies in heavy nuclei. 

The results of TD RMF calculations for ^'^^Pb are displayed in Fig. 4: time-dependent 
monopole moments = ^($(t)|r^|$(t)) and the corresponding Fourier power spectra 

for the nonlinear effective interactions. As one would expect for a heavy nucleus, there 
is very little spectral fragmentation and a single mode dominates, at least for NLl and 
NL3. The experimental excitation energy 13.7 ± 0.3 McV is very close to the frequency 
of oscillations obtained with the NL3 parameter set: 14.1 McV. The calculated excitation 
energy for the NLl parameter set (iCnm = 211.7 MeV), is approximately 1 MeV lower than 
the average experimental value. For the linear effective forces HS and LI the oscillations 
are more anharmonic, and the monopole strength is located well above the experimental 
GMR energy. The effective interactions NLl and NL3 seem to produce GMR excitation 
energies which are quite close to the experimental values. For these two parameter sets we 
have calculated the isoscalar giant monopole resonances in a number of doubly closed-shell 
nuclei: ^°Ca, ^^Ni, ioo,ii4,i32gj^^ ^^'^'^'^Zy, ^^^Gd. The results are shown in Fig. 5. The energies 
of giant monopole states are determined from the Fourier spectra of the time-dependent 
monopole moments, and are displayed as function of the mass number. The NLl excitation 
energies arc systematically lower, but otherwise the two effective interactions produce very 
similar dependence on the mass number. The empirical curve 

included in the figure, and it follows very closely the excitation energies calculated with the 
NL3 parameter set. Similar results are obtained from constrained RMF calculations. Both 

methods indicate that, in the framework of relativistic mean field theory, the nuclear matter 
compression modulus ~ 250 — 270 MeV is in reasonable agreement with the available 
data on spherical nuclei. This value is approximately 20% larger than the values deduced 
from recent non-relativistic density dependent Hartree-Fock calculations with Skyrme or 
Gogny forces 

2.3 Regular and chaotic dynamics of collective vibrations 

The atomic nucleus has been used as a laboratory, both experimentally and theoretically, 
for investigating the transition from order to chaos in quantum dynamical systems. Most of 
these studies have concentrated on two major aspects: (i) generic signatures of chaos in local 
fiuctuations and correlations of nuclear level distributions, and (ii) chaos in microscopic and 
collective dynamics of realistic many-body systems. Regular and chaotic dynamics in giant 
nuclear oscillations has been the subject of a number of studies. What has emerged as a 
very interesting result is that an undamped collective mode may coexists with chaotic single- 
particle motion. It appears that the slowly vibrating self-consistent mean field created by 
the same nucleons averages out the random components in their motion. In all investigations 
the motion of only one type of particles has been considered. That is, only the dynamics of 
isoscalar collective modes. 

We have studied the difference in the dynamics of isoscalar and isovector collective modes. 



6 



In particular, we consider isoscalar and isovector monopole oscillations in spherical nuclei, 
but analogous considerations apply to higher multipolarities. In Fig. 6 results are shown of 
time-dependent relativistic mean-field calculations for isoscalar and isovector oscillations in 
^'^^Pb. The experimental isoscalar GMR energy in ^°^Pb is 13.7±0.3 MeV, and the excitation 
energy of the isovector mode is 26 ± 3 MeV. Calculations have been performed for the NLl 
effective interaction. In the isoscalar case both proton and neutron densities are radially 
expanded, while for the isovector mode the proton density is initially compressed by the same 
amount. Therefore, in both cases we follow the time evolution of the same system, just the 
initial conditions are different. In Fig. 6a we plot the time history of the isoscalar monopole 
moment (r^(t)), and in Fig. 6b the corresponding isovector moment < r'^{t) > — < r'^{t) > is 
displayed. The isoscalar mode displays regular undamped oscillations, while for the isovector 
mode we observe strongly damped anharmonic oscillations. On the right hand panels we 
plot the corresponding Fourier power spectra. The Fourier spectrum of the isovector mode 
is strongly fragmented. However, the main peaks are found in the energy region 25 — 30 
MeV, in agreement with the experimental data. For the isoscalar mode, the time history of 
the monopole moment and the Fourier spectrum show that the oscillations of the collective 
coordinate are regular. On the other hand, the appearance of a broad spectrum of frequencies 
seems to indicate that the isovector oscillations are chaotic. 

A diagnosis of chaotic vibrations would imply that one has a clear definition of such motion. 
For a quantum system, however, the concept of chaos, especially in time-dependent prob- 
lems, is not well defined. And although our description of nuclear vibrations is semi-classical, 
quantum effects like the Pauli principle are present in the initial conditions and during the 
dynamical evolution. There exists a number of tests that can help to identify chaotic oscilla- 
tions in physical systems, and some of them can be applied in the present consideration. In 
Figs. 7 - 9 we display some additional qualitative measures which can be used to characterize 
the response of our nonlinear system. In Fig. 7 we have constructed the two-dimensional 
time-delayed pseudo-phase space for the isoscalar (a) and isovector (b) oscillations, shown 
in Fig. 6. Since information is available on the time evolution of just one variable, the 
collective coordinate, one plots the signal versus itself, but delayed or advanced by a fixed 
time constant [< r^(t) >, < r^(t-|-r) >]. The phase space trajectories for the isoscalar mode 
are closed ellipses, indicating regular oscillations. For the isovector oscillations on the other 
hand, the trajectories are completely chaotic. The strong damping results from one-body 
processes: (i) escape of nucleons into the continuum states and (ii) coUisions of the nucleons 
with the moving wall of the nuclear potential generated by the self-consistent mean fields. 
In Fig. 8 we display the corresponding Poincare sections constructed from three-dimensional 
time-delayed pseudo-phase space. The Poincare map for the isoscalar mode consists of two 
sets of closely located points, and therefore confirms regular oscillations. For the isovector 
oscillations the Poincare map appears as a cloud of unorganized points in the phase plane. 
Such a map indicates stochastic motion. Another measure that is related to the Fourier 
transform is the autocorrelation function 



T 




(11) 
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When the signal is chaotic, information about its past origins is lost. This means that 
A{t) ^ as r ^ oo, or the signal is only correlated with its recent past. The autocorrelation 
functions for isoscalar and isovector oscillations are shown in Fig. 9. For the isovector mode 
A{t) displays a rapid decrease, and the envelope appears as an irregular waveform. The 
oscillations of the collective coordinate can be characterized as regular for the isoscalar 
mode, and they become chaotic when initial conditions correspond to the isovector mode. In 
Ref. we have shown how a regular collective mode can coexist with chaotic single-particle 
dynamics, a result which confirms the conclusions of a number of studies. However, we have 
also shown that this is the case only for isoscalar modes, that is, only if one considers the 
motion of a single type of particles. When protons and neutrons move out of phase, as 
it happens for isovector modes, the resulting dynamics of the collective coordinate exhibit 
chaotic behavior. This is explained by the fact that protons and neutrons effectively move in 
two self-consistent potentials that oscillate out of phase. For example, when neutrons move 
inward, they scatter on the potential wall with positive curvature that is created by protons 
moving outward. This will lead to pseudo-random motion of the nucleons and dissipation of 
collective oscillations. 
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Figure Captions 

FIG. 1. Time-dependent quadrupole moment and the corresponding Fourier spectrum for 
^°Ca. 

FIG. 2. Fourier spectra of the time-dependent baryon density for various values of the 
coordinate on the axis along which the initial densities are deformed. 

FIG. 3. Mean-field action as function of the excitation energy of the nucleons E* —< 
^{t)\HDit)\^{t) > - < ^Gs\HDiO)\^GS > (a), and of the initial deformation P (b). 

FIG. 4. Time-dependent isoscalar monopole moments < > (t) and the corresponding 
Fourier power spectra for ^°^Pb. The parameter sets are NLl, NL3, NL-SH and NL2. 

FIG. 5. Excitation energies of isoscalar giant monopole resonances in spherical nuclei as 
function of the mass number. The effective interactions are: NLl (squares) and NL3 
(circles). The solid curve corresponds to the empirical relation 80 A~^^^ MeV. 

FIG. 6. Results of time-dependent relativistic mean-field calculations for isoscalar and 
isovector oscillations in ^°^Pb. 

FIG. 7. Pseudo-phase space for isoscalar (a) and isovector (b) monopole oscillations in 

208pb. 

FIG. 8. Poincare sections for isoscalar (a) and isovector (b) monopole oscillations in ^°^Pb. 
FIG. 9. Autocorrelation functions for isoscalar and isovector oscillations in ^°^Pb. 
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Abstract 

The relativistic mean-field theory provides a framework in which the nuclear many- 
body problem is described as a self-consistent system of nucleons and mesons. In 
the mean-field approximation, the self-consistent time evolution of the nuclear system 
describes the dynamics of collective motion: double giant resonances, nuclear com- 
pressibility from monopole resonances, regular and chaotic dynamics of isoscalar and 
isovector collective vibrations. 



1 Introduction and outline of the model 

Relativistic mean-field (RMF) models have been successfully applied in calculations of nu- 
clear matter and properties of finite nuclei throughout the periodic table. In the self- 
consistent mean-field approximation, detailed calculations have been performed for a variety 
of nuclear structure phenomena [1]. In the present work we review the applications of RMF 
to the dynamics of collective vibrations in spherical nuclei. In relativistic quantum hadro- 
dynamics [2], the nucleus is described as a system of Dirac nucleons which interact through 
the exchange of virtual mesons and photons. The Lagrangian density of the model is 

£ = ^(«^.a-m)V' + ]^{daf-U{a) 

- gai^o-tjj - g^t/j-f ■ ujtjj - gpt/j-f ■ prt/j - et/j-f ■ A^^—^-^t/j . (1) 

The Dirac spinor t/j denotes the nucleon with mass m. m(j, m^^, and nip are the masses of the 
cr-meson, the cu-meson, and the /j-meson, and g^j^ g^^^ and gp are the corresponding coupling 
constants for the mesons to the nucleon. U (cr) denotes the nonlinear a self-interaction 

U{a) = ^m'y + ^g2a^ + jg3<T\ (2) 
^Alexander von Humboldt Fellow, on leave of absence from University of Zagreb, Croatia 
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and 0^^ and F^"" are field tensors [2]. 

The coupled equations of motion are derived from the Lagrangian density (1). The Dirac 
equation for the nucleons: 



OL 



(1 -Ts) , 
+5'c^^o + gpTpo + e Ao 



2 



and the Klein- Gordon equations for the mesons: 



-A + ml ) Lo 
- A + m„ , 



p 

a1 a 



9pJ p 



P 



p 



(3) 



(4) 
(5) 
(6) 
(7) 



In the relativistic mean-field approximation, the nucleons described by single-particle spinors 
t/ji [i = 1,2,..., A) are assumed to form the A-particle Slater determinant |$), and to 
move independently in the classical meson fields. The sources of the fields, i.e. den- 
sities and currents, are calculated in the no-sea approximation [3]: the scalar density: 



Ps = 

bar yon current: j 



J2i=i i^ii^i^ the isoscalar baryon current: j' 



1^8=1 V'j'T^V'i? the isovector 



nP 

J era 



= J2i=i ''Pn^^i^ij the electromagnetic current for the photon-field: 
J2t=i V'«'7^''~^V'«- The summation is over all occupied states in the Slater determinant 
1$). Negative-energy states do not contribute to the densities in the no-sea approximation of 
the stationary solutions. It is assumed that nucleon single-particle states do not mix isospin. 
The ground state of a nucleus is described by the stationary self-consistent solution of the 
coupled system of equations (3)-(7), for a given number of nucleons and a set of coupling 
constants and masses. The solution for the ground state specifies part of the initial conditions 
for the time-dependent problem. The dynamics of nuclear collective motion is analyzed in 
the framework of time-dependent relativistic mean-field model, which represents a relativis- 
tic generalization of the time-dependent Hartree-Fock approach. For a given set of initial 
conditions, i.e. initial values for the densities and currents, nuclear dynamics is described by 
the simultaneous evolution of A single- particle Dirac spinors in the time-dependent mean 
fields. Frequencies of eigenmodes are determined from a Fourier analysis of dynamical quan- 
tities. In this microscopic model, self-consistent time-dependent mean-field calculations are 
performed for multipole excitations. An advantage of the time-dependent approach is that 
no assumption about the nature of a particular mode of vibrations has to be made. Retar- 
dation effects for the meson fields are not included in the model, i.e. the time derivatives 
in the equations of motions for the meson fields are neglected. This is justified by the large 
masses in the meson propagators causing a short range of the corresponding meson exchange 
forces. Negative energy contributions are included implicitly in the time-dependent calcu- 
lation, since the Dirac equation is solved at each step in time for a different basis set [3]. 
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Negative energy components with respect to the original ground-state basis are taken into 
account automatically, even if at each time step the no-sea approximation is applied. 
The description of nuclear dynamics as a time-dependent initial- value problem is intrinsi- 
cally semi-classical, since there is no systematic procedure to derive the initial conditions 
that characterize the motion of a specific mode of the nuclear system. The theory can be 
quantized by the requirement that there exist time-periodic solutions of the equations of 
motion, which give integer multiples of Planck's constant for the classical action along one 
period [4]. For giant resonances the time-dependence of collective dynamical quantities is 
actually not periodic, since generally giant resonances are not stationary states of the mean- 
field Hamiltonian. The coupling of the mean-field to the particle continuum allows for the 
decay of giant resonances by direct escape of particles. In the limit of small amplitude os- 
cillations, however, the energy obtained from the frequency of the oscillation coincides with 
the excitation energy of the collective state. In Refs. [3, 4, 5] we have shown that the model 
reproduces experimental data on giant resonances in spherical nuclei. 

2 Dynamics of collective vibrations 
2.1 Double Giant Resonances 

The physics of giant resonances in nuclei has been the subject of extensive experimental 
and theoretical studies for many years. However, only recently substantial evidence has 
been reported for two-phonon states built with giant resonances. Resonant structures that 
were observed in heavy-ion inelastic scattering, have been interpreted as possible multiple 
excitations of the giant quadrupole resonance. Double giant dipole resonances have also been 
discovered in the neutron and in the 7-spectra of nuclei that have been Coulomb excited in 
relativistic heavy-ion collisions. In Ref. [4] we have performed time-dependent relativistic 
mean-field calculations and found evidence for modes which can be interpreted as double 
resonances, and which in a quantized theory correspond to two-phonon states. 
As an example of double resonances in light nuclei, we consider the double isoscalar giant 
quadrupole resonance in "^^Ca. The experimental spectrum exhibits a prominent structure, 
centered at 34 it 2 MeV excitation energy, with a width of 9 it 2 MeV. It is interpreted as 
the two-phonon state of the single isoscalar GQR at 17.5 MeV. Using the NL-SH parameter 
set for the effective Lagrangian, we have studied isoscalar quadrupole oscillations in "^^Ca. 
The quadrupole mode of oscillations is excited by deforming the spherical solution for the 
ground-state. For a specific initial deformation, we follow the time-evolution of the collective 
variable, the quadrupole moment 

q2o{t) = {m\Q2o\m) = {m\r'Y2o\m) (s) 

The time-dependent quadrupole moment shown in Fig. 1 corresponds to an initial axial 
deformation of the baryon density /3 = 0.38. The resulting Fourier spectrum displays a 
strong peak at 18.5 MeV, in reasonable agreement with the experimental data. 
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We find evidence for excitation of a liiglier mode in tlie oscillations of the baryon density. 
The wave function of the nuclear system is a Slater determinant at all times, and therefore 
can be expanded in the basis of the ground state |$o) 

= |$o) I^o) + E ^™(0^™v(t)a+a+,a,,a, |$o) + ... (9) 

mi mim'i' 

If the total wave function contains collective 2p — 2h components, they will be observed in 
the Fourier spectrum of the time-dependent baryon density Pb{^^ t) = ZliLi ''Pti^^ OV'^l'"? 0- 
Because of axial symmetry and the isoscalar nature of the excitation, it is sufficient to 
consider oscillations in time of the baryon density on the positive 2;-axis. In Fig. 2 we 
display the Fourier transforms of the time-dependent baryon density for various values of 
the coordinate z. Two peaks are clearly observed. The first one at 18.5 MeV corresponds to 
the isoscalar quadrupole resonance. It gradually increases from the center toward the surface 
of the nucleus. If we plot the values of the Fourier transforms at 18.5 MeV as a function of 
the resulting curve corresponds to the transition density. The transition density for the first 
peak is typical for isoscalar quadrupole resonances. The second peak is at 37 MeV, twice 
the energy of the GQR. It has a maximum in the center of the nucleus, at first decreases 
with z, but then appears again on the surface. Compared to the GQR transition density, 
the curve for the 37 MeV peak displays an additional node. 

Periodic solutions of the time-dependent Dirac equations can be used to construct the ener- 
gies of the many-body system. The energy spectrum can be obtained from a semi-classical 
quantization procedure. One finds periodic solutions such that the mean-field action along 
a periodic orbit 



I = E dt 



to 



d 

< 'ijj,{t)\ifi—\'ijj,{t) > - e, 



(10) 
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is equal to an integer multiple of the Planck constant / = nh. In Eq. (10) T is the period 
of oscillations, \4^i{t) > denotes time-dependent single-nucleon Dirac spinors, are the 
corresponding single-nucleon energies in the unperturbed ground-state, and the summation 
runs over the occupied states. Because of the coupling to the continuum in the mean-field 
description, giant resonances are not stationary states of the Hamiltonian. Consequently, 
a non-periodic dependence on time is obtained for dynamical quantities. If the damping is 
very strong, the giant resonance is not periodic even on the average, and the quantization 
condition cannot be applied. However, if the motion is nearly periodic, i.e. the damping is 
relatively weak, the quantization procedure can still be used to calculate the energies, and 
the effect of damping can be taken into account approximately. 
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FIG. 2. Fourier spectra of the time-dependent baryon density for various values of the coordinate 
z, on the axis along which the initial densities are deformed. 



In Fig. 3 we display the mean-field action as function of the excitation energy of the nucleons 
and of the initial deformation /3. For the values of /3 indicated by dots we have integrated 
the coupled system of Dirac and Klein-Gordon equations and calculated the action integral. 
The mean-field action is a quadratic function of the initial deformation /3, and an almost 
perfect linear function of the excitation energy. Only above 60 MeV a slight deviation from a 
pure linear dependence is observed. The values of < £"* > for which the action is an integer 
multiple of the Planck constant are: I = 1 h for 18.5 MeV, and I = 2 h for 37.1 MeV. 
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FIG. 3. Mean-field action as function of the excitation energy of the nucleons 

E* =< ^{t)\HD{t)\^{t) > - < ^Gs\HDiO)\^GS > (a), and of the initial deformation (3 (b). 

A double giant dipole resonance has been observed in relativistic Coulomb excitation of 
208pb. The single GDR is found at 13.3 ± 0.1 MeV with a width of 4.1 ± 0.1 MeV. The 
sum energy of coincident photon pairs displays a broad structure at 25.6 it 0.9 MeV with a 
width of 5.8 lb 1.1 MeV. It is interpreted as the double GDR. In order to excite isovector 
dipole motion we define the initial conditions: at t = (in the center of mass system) all 
protons start moving in the -\-z direction with velocity u^, and all neutrons start moving 
in the —z direction with velocity = j^v^^. For the NL-SH parameter set, the Fourier 
spectrum of the time-dependent dipole moment displays a strong peak at 12.9 MeV excitation 
energy, in good agreement with experimental data. We have found that the mean-field 
action is an integer multiple of the Planck constant: I = 1 h for < E* >= 12.9 MeV, and 
I = 2 h for < E* >= 25.9 MeV. Therefore, the energy of the one-phonon state, calculated 
from the mean-field action, coincides with the resonant energy of the mean peak in the 
Fourier spectrum, and the two-phonon state at 25.9 MeV is in excellent agreement with the 
experimental value for the excitation energy of the double GDR. 

2.2 Monopole Giant Resonances and nuclear compressibility 

The study of isoscalar monopole resonances in nuclei provides an important source of in- 
formation on the nuclear matter compression modulus Knm- This quantity is crucial in 
the description of properties of nuclei, supernovae explosions, neutron stars, and heavy-ion 
collisions. In principle the value of Anm can be extracted from experimental energies of 
isoscalar monopole vibrations in nuclei (giant monopole resonances GMR). However, the 
complete experimental data set on isoscalar GMR does not limit the range of Knm to better 
than 200 — 350 MeV. Microscopic calculations of GMR excitation energies might provide 
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a more reliable approach to the determination of the nuclear matter compression modulus. 
Modern non-relativistic Hartree-Fock plus RPA calculations, using both Skyrme and Gogny 
effective interactions, indicate that the value of Knm should be in the range 2f0-220 MeV. 
In relativistic calculations on the other hand, both time-dependent and constrained RMF 
results indicate that empirical GMR energies are best reproduced by an effective force with 
A^nm ~ 250 - 270 MeV. 
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FIG. 4. Time-dependent isoscalar monopole moments < > (t) and the corresponding Fourier 
power spectra for ^°®Pb. The parameter sets are NLl, NL3, NL-SH and NL2. 



In Ref. [5] we have performed time-dependent and constrained RMF calculations for 
monopole giant resonances of a number of spherical closed shell nuclei, from ^^0 to the 
heavy nucleus ^°®Pb. For the effective Lagrangian we have used six parameter sets, which 
differ mostly by their prediction for A^nm, but otherwise reproduce reasonably well experi- 
mental data on nuclear properties. The idea is to restrict the possible values of the nuclear 
matter compression modulus, on the basis of the excitation energies of giant monopole states 
calculated with different effective interactions. In addition to the four non-linear sets NLI, 
NL3, NL-SH and NL2, we have also included two older linear parametrizations, HS and LI. 
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The sets NLl, NL-SH and NL2 have been extensively used in the description of properties 
of finite nuclei [1]. In order to bridge the gap between NLl (A^nm = 211.7 MeV), and NL-SH 
(-Anm = 355.0 MeV), we have also included a new effective interaction NL3 [6] (A^nm = 271.8 
MeV). This new parameter set provides an excellent description not only for the properties of 
stable nuclei, but also for those far from the valley of beta stability. From the energy spectra 
and transition densities calculated with these effective forces, it has been possible to study 
the connection between the incompressibility of nuclear matter and the breathing mode en- 
ergy of spherical nuclei. For the isoscalar mode we have found an almost linear relation 
between the excitation energy of the monopole resonance and the nuclear matter compres- 
sion modulus. For the determination of Knm especially relevant are microscopic calculations 
of GMR excitation energies in heavy nuclei. The results of TD RMF calculations for ^°®Pb 
are displayed in Fig. 4: time-dependent monopole moments (r^(t)) = ^($(t)|r^|$(t)) 
and the corresponding Fourier power spectra for the nonlinear effective interactions. As one 
would expect for a heavy nucleus, there is very little spectral fragmentation and a single 
mode dominates, at least for NLl and NL3. The experimental excitation energy 13.7 it 0.3 
MeV is very close to the frequency of oscillations obtained with the NL3 parameter set: 14.1 
MeV. The calculated excitation energy for the NLl parameter set (A^nm = 211.7 MeV), is 
approximately 1 MeV lower than the average experimental value. For the linear effective 
forces HS and LI the oscillations are more anharmonic, and the monopole strength is located 
well above the experimental GMR energy. 

20 1- 1 ■ 1 ' < 1 




12 1 ■ 1 • u_ 
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FIG. 5. Excitation energies of isoscalar giant monopole resonances in spherical nuclei as function 
of the mass number. The effective interactions are: NLl (squares) and NL3 (circles). The solid 
curve corresponds to the empirical relation ~ 80 A~^^'^ MeV. 

The effective interactions NLl and NL3 seem to produce GMR excitation energies which 
are quite close to the experimental values. For these two parameter sets we have calculated 
the isoscalar giant monopole resonances in a number of doubly closed-shell nuclei: "^^Ca, 
s^Ni, i°°'ii4'i32gj^^ 90'i222r^ i46Qj_ rj,^^ rcsults are shown in Fig. 5. The energies of giant 
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monopole states are determined from the Fourier spectra of the time-dependent monopole 
moments, and are displayed as function of the mass number. The NLl excitation energies 
are systematically lower, but otherwise the two effective interactions produce very similar 
dependence on the mass number. The empirical curve 80 jvieV is also included 

in the figure, and it follows very closely the excitation energies calculated with the NL3 pa- 
rameter set. Similar results are obtained from constrained RMF calculations. Both methods 
indicate that, in the framework of relativistic mean field theory, the nuclear matter com- 
pression modulus Knm ~ 250 — 270 MeV is in reasonable agreement with the available data 
on spherical nuclei. This value is approximately 20% larger than the values deduced from 
recent non-relativistic density dependent Hartree-Fock calculations with Skyrme or Gogny 
forces 

2.3 Regular and chaotic dynamics of collective vibrations 

The atomic nucleus has been used as a laboratory, both experimentally and theoretically, 
for investigating the transition from order to chaos in quantum dynamical systems. Most of 
these studies have concentrated on two major aspects: (i) generic signatures of chaos in local 
fiuctuations and correlations of nuclear level distributions, and (ii) chaos in microscopic and 
collective dynamics of realistic many-body systems. Regular and chaotic dynamics in giant 
nuclear oscillations has been the subject of a number of studies. What has emerged as a 
very interesting result is that an undamped collective mode may coexists with chaotic single- 
particle motion. It appears that the slowly vibrating self-consistent mean field created by 
the same nucleons averages out the random components in their motion. In all investigations 
the motion of only one type of particles has been considered. That is, only the dynamics of 
isoscalar collective modes. 

We have studied the difference in the dynamics of isoscalar and isovector collective modes. 
In particular, we consider isoscalar and isovector monopole oscillations in spherical nuclei, 
but analogous considerations apply to higher multipolarities. In Fig. 6 results are shown of 
time-dependent relativistic mean-field calculations for isoscalar and isovector oscillations in 
^°®Pb. The experimental isoscalar GMR energy in ^°®Pb is I3.7ib0.3 MeV, and the excitation 
energy of the isovector mode is 26 it 3 MeV. Calculations have been performed for the NLI 
effective interaction. In the isoscalar case both proton and neutron densities are radially 
expanded, while for the isovector mode the proton density is initially compressed by the same 
amount. Therefore, in both cases we follow the time evolution of the same system, just the 
initial conditions are different. In Fig. 6 we plot the time history of the isoscalar monopole 
moment (r^(t)), and in Fig. lb the corresponding isovector moment < r^{t) > — < r^{t) > is 
displayed. The isoscalar mode displays regular undamped oscillations, while for the isovector 
mode we observe strongly damped anharmonic oscillations. On the right hand panels we 
plot the corresponding Fourier power spectra. The Fourier spectrum of the isovector mode 
is strongly fragmented. However, the main peaks are found in the energy region 25 — 30 
MeV, in agreement with the experimental data. For the isoscalar mode, the time history of 
the monopole moment and the Fourier spectrum show that the oscillations of the collective 
coordinate are regular. On the other hand, the appearance of a broad spectrum of frequencies 
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seems to indicate that the isovector oscillations are chaotic. 



FIG. 6. 
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A diagnosis of chaotic vibrations would imply that one has a clear definition of such motion. 
For a quantum system, however, the concept of chaos, especially in time-dependent prob- 
lems, is not well defined. And although our description of nuclear vibrations is semi-classical, 
quantum effects like the Pauli principle are present in the initial conditions and during the 
dynamical evolution. There exists a number of tests that can help to identify chaotic oscilla- 
tions in physical systems, and some of them can be applied in the present consideration. In 
Figs. 7-9 we display some additional qualitative measures which can be used to characterize 
the response of our nonlinear system. In Fig. 7 we have constructed the two-dimensional 
time-delayed pseudo-phase space for the isoscalar (a) and isovector (b) oscillations, shown 
in Fig. 6. Since information is available on the time evolution of just one variable, the 
collective coordinate, one plots the signal versus itself, but delayed or advanced by a fixed 
time constant [< r^(t) >, < r^(t-|-T) >]. The phase space trajectories for the isoscalar mode 
are closed ellipses, indicating regular oscillations. For the isovector oscillations on the other 
hand, the trajectories are completely chaotic. The strong damping results from one-body 
processes: (i) escape of nucleons into the continuum states and (ii) collisions of the nucleons 
with the moving wall of the nuclear potential generated by the self-consistent mean fields. 
In Fig. 8 we display the corresponding Poincare sections constructed from three-dimensional 
time-delayed pseudo-phase space. The Poincare map for the isoscalar mode consists of two 
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sets of closely located points, and therefore confirms regular oscillations. For the isovector 
oscillations the Poincare map appears as a cloud of unorganized points in the phase plane. 
Such a map indicates stochastic motion. 
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FIG. 7. Pseudo-phase space for isoscalar (a) and isovector (b) monopole oscillations in Pb 
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FIG. 8. Poincare sections for isoscalar (a) and isovector (b) monopole oscillations in Pb. 
Another measure that is related to the Fourier transform is the autocorrelation function 



A{t) = lim / < r^{t) >< r^{t + r) > dt 





'11] 



When the signal is chaotic, information about its past origins is lost. This means that 
A(t) — )• as t — )• oo, or the signal is only correlated with its recent past. The autocorrelation 
functions for isoscalar and isovector oscillations are shown in Fig. 9. For the isovector mode 
A(t) displays a rapid decrease, and the envelope appears as an irregular waveform. 
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The oscillations of the collective coordinate can be characterized as regular for the isoscalar 
mode, and they become chaotic when initial conditions correspond to the isovector mode. In 
Ref. [7] we have shown how a regular collective mode can coexist with chaotic single-particle 
dynamics, a result which confirms the conclusions of a number of studies. However, we have 
also shown that this is the case only for isoscalar modes, that is, only if one considers the 
motion of a single type of particles. When protons and neutrons move out of phase, as 
it happens for isovector modes, the resulting dynamics of the collective coordinate exhibit 
chaotic behavior. This is explained by the fact that protons and neutrons effectively move in 
two self-consistent potentials that oscillate out of phase. For example, when neutrons move 
inward, they scatter on the potential wall with positive curvature that is created by protons 
moving outward. This will lead to pseudo-random motion of the nucleons and dissipation of 
collective oscillations. 
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